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Abstract 

Individual taste differences were first reported in the first half of the 20th 
century, but the primary reasons for these differences have remained uncertain. Much of 
the taste variation among different mammalian species can be explained by 
pseudogenization of taste receptors. In this study, by analyzing 14 ethnically diverse 
populations, we investigated whether the most recent disruptions of taste receptor genes 
segregate with their intact forms. Our results revealed an unprecedented prevalence of 
segregating loss-of-function (LoF) taste receptor variants, identifying one of the most 
pronounced cases of functional population diversity in the human genome. LoF variant 
frequency (2.10%) was considerably higher than the overall mutation rate (0.16%), and 
many humans harbored varying numbers of critical mutations. In particular, molecular 
evolutionary rates of sour (14.7%) and bitter receptors (1.8%) were far higher in 
humans than those of sweet (0.02%), salty (0.05%), and umami (0.17%) receptors 
compared with other carnivorous mammals although not all of the taste receptors genes 
were identified. Many LoF variants are population-specific, some of which arose even 
after the population differentiation, but not before divergence of the modern and archaic 
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(Neanderthal and Denisovan) human. Based on these findings, we conclude that modern 
humans might have been losing their taste receptor genes because of high-frequency 
LoF taste receptor variants. Finally I actually demonstrated the genetic testing of taste 
receptors from personal exome sequence. 

Introduction 

Taste is a primitive sense that helps recognize and distinguish key dietary 
components (1-3). Taste perception affect nutrient intake and appetite, and thus taste 
decrements can lead to food poisoning and over-exposure to hazardous chemicals 
(1,4). The sensation of taste is classified into five prototypical categories: sweet, bitter, 
sour, salty, and umami (1-3,5-7). In addition recent compelling evidence raises the 
possibility of an additional sixth and seventh taste modality devoted to the perception of 
calcium (8,9) and lipids (10,11). We all perceive these tastants through sensory organs 
called taste buds (2,3,5-7). Humans have approximately 10,000 taste buds, primarily 
located on the tongue with a few on the palate and throat (1). 

The anatomical units of taste detection are taste receptor cells (TRCs) (1-3). 

TRCs are distributed among different lingual papillae and palate epithelia, and are 
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involved in recognizing various taste stimuli (1-3). The sensation of taste is produced 
when stimuli in the mouth react with specific taste receptors (1-3,5-7). Taste coding at 
the peripheral cells may rely on straightforward labeled lines (that is, sweet, bitter, sour, 
salty, and umami signals) to transform tastant quality into neural signals through the 
glossopharyngeal and facial nerve (1,3). Taste in conjunction with olfaction is perceived 
as flavors or other substances (1-3). 

During the last 10 years, the discovery and characterization of mammalian taste 
receptors has progressed tremendously (1-3,6,7). The recent identification of cells and 
receptors that mediate the prototypical tastes has generated powerful molecular tools 
that can be used to devise tests to establish the basis of taste coding at the periphery 
(12-18); tastant-induced activity were measured using nerves innervating the buds in 
various gene-knockout (KO) mice (2,15) or receptor-expressing cells (16-18), and 
diphtheria toxin A fragment (DTA)-mediated cell ablation studies were performed to 
uncover the specific taste stimuli in vivo (19). 

"Why do some people have different taste recognition thresholds compared 
with others?" "Do they have different receptor densities within their body?" These have 
been controversial issues since the early 20th century and the reasons for individual 

4 



Downloaded from http://biorxiv.org/on September 18, 2014 



taste differences remain uncertain (20-26). Recent studies have raised the possibility 
that most genotype-phenotype association studies are based on rare gene disruptions 
arising during mammalian evolution (27,28). In fact, some taste receptors also seem to 
be affected by mutations during mammalian evolution (6,7,12,13,29-31). For example, 
the hT2R43 W35 allele makes people most sensitive to the bitterness of the natural plant 
toxins, aloin and aristolochic acid, and an artificial sweetener, saccharin, compared with 
S35 allele (31). In addition, a closely related gene's HT2R44 W35 allele is also sensitive 
to saccharin (31). People who do not possess these allele do not taste the bitterness of 
these compounds at low micromolar concentrations (31). Another well-known example 
is KTAS2R38. African populations have higher levels of genetic, geographic and 
phenotypic diversity at the TAS2R38 locus (32). TAS2R38 P49A, A262V and V296I 
have clear effects on the phenylthiocarbamide (PTC) sensitivity in humans (33-35). 
These polymorphisms illustrate the influence of recent genetic variation on a common 
trait. However none of the genetic variants explain the whole spectrum of 
interindividual differences of taste sensitivity. 

We postulated that individual taste differences were predominantly caused by 

loss-of-function (LoF) variants of taste receptor genes. In this study, we investigated 

whether the most recent of these disruptions may still segregate with the intact forms by 
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analyzing 14 ethnically diverse populations. The recently sequenced completed 
genomes for many modern humans provided the opportunity for reconstruction of the 
currently known LoF variants (37,38). Our results revealed an unprecedented 
prevalence of segregating LoF variants, one of the most pronounced cases of functional 
population diversity in the human genome. 

Results 

Profiles of human taste receptor gene variations reveal striking individuality 

The representative taste receptors are composed of more than 50 coding 

regions distributed in clusters over most chromosomes in mammals (5). In mice, taste 

receptor pseudogenes comprise 15% of this gene range; in humans, a fraction roughly 

two times larger appears to be inactivated (5). Of more than 50 taste receptor genes, 

more than half appear to have been pseudogenized by mutations during mammalian 

evolution (5). Extreme diminution of the functional taste receptor repertoire was a 

relatively recent genomic process and is probably still ongoing (5). Therefore, we 

conjectured that a substantial fraction of modern human taste receptors may segregate 
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between intact and pseudogene forms. Isolated cases of segregating mutations in taste 
receptors have been reported (6,7,12,13,29-31). 

We focused on three kinds of LoF variants expected to correlate with complete 

loss of function of the affected transcripts: one stop codon-introducing (nonsense) or 

two (5' and 3') splice site disrupting single-nucleotide substitutions. We searched 1,000 

genomes (37) and NHLBI databases (38) for variations with the potential to affect the 

protein integrity of taste receptor genes (fig. 1 and table SI). Target receptor genes are 

as follows: TAS1R1, TAS1R2, TAS1R3, TAS2R1, TAS2R3, TAS2R4, TAS2R5, TAS2R7, 

TAS2R8, TAS2R9, TAS2R10, TAS2R13, TAS2R14, TAS2R16, TAS2R19, TAS2R20, 

TAS2R30, TAS2R31, TAS2R38, TAS2R39, TAS2R40, TAS2R41, TAS2R42, TAS2R43, 

TAS2R45, TAS2R46, TAS2R50, TAS2R60, PKD1L3, PKD2L1, CD36, ENaCa, ENaCd, 

HCN1, HCN4, TRPM5, and CALHM1 (table SI). I did not analyze the putative taste 

receptor genes and perform structural similarity searches. Taste receptor are different 

from olfactory receptors (ORs) in that taste receptors are not composed of single gene 

family (1-3). Our data matrix included approximately 8,000 individuals from six ethnic 

origins (table S2). In this largest realm, we extrapolated that the number of major 

segregating taste receptor LoF variants in the entire human genome had to be at least 24, 

of which 18 were expected to have a major allele frequency of >0.05% (Figure 2A). The 
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average LoF variant frequency among known taste receptor genes (2.10%; TAS2R (18) 
= 1.76%; PKD-like (PKD 1 L3+PKD2L1 (2,3) = 14.7%) was extremely higher than the 
overall mutation rate in the human genome (approximately 0.16% (calculated from LoF 
variants caused by nucleotide substitution; this proportion also includes both LoF 
frequency of taste and smell receptor genes) (38)) (fig. 2B). Except for complete 
pseudogenes, mutation rates of taste receptor genes exceeded those of olfactory receptor 
(OR) genes (approximately 1.41%; recalculated from 367 functional ORs because more 
than half genes analyzed in a past study (39) have already been proved to be complete 
pseudogenes; If we calculate the mutation frequency of pseudogenes, their MAF is 
regarded as 100% and mutation frequency become quite high). 

In theory, a combination of these genes could give rise to an enormous series of 
individual differences in taste perception; each examined individual had a unique 
genotypic pattern. The average number of LoF sites in TAS1R (2,3), TAS2R (18) and 
PKD-like (PKD 1 L3+PKD2L1 (2,3)) receptors per individual was 0.00371, 0.422 and 
0.294, respectively. Coupled with TAS and PKD-like genes, CD36 (fat receptor 
(10,11)), and ENaCd (SCNN1D) (sodium channel (40-43))) were frequently lost in 
some populations (fig. 2A, table S2 and S3). In contrast, we confirmed that the other 
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taste buds-specific or -enriched genes previously reported (44,45) had been rarely lost 
(forty five genes analyzed are described in supplementary note). 

No novel taste receptor genes were acquired after the divergence of 
chimpanzees and humans. Using 1000 genome sequences, we confirmed that the three 
pseudogenized bitter receptors TAS2R2P, TAS2R62P, and TAS2R64P (46) have never 
been functionally restored by gain-of-function mutations in any human population. Our 
survey of LoF variants in taste receptor genes pointed out the possibility that modern 
humans might have continued to genetically lose the repertoire of receptor genes after 
species differentiation. 

Population-specific LoF variants in taste receptor genes 

Because of physical barriers to migration, ethnic populations rarely interbreed as 
convenient in theoretical random models. Consequently populations from different 
ethnic groups often have different genetic backgrounds and therefore different 
frequencies of genetic polymorphisms (47). To examine the genetic distance of taste 
receptor genes among populations, we first investigated whether frequencies and 
distributions of mutations had regional and racial biases by constructing a map for 

9 



Downloaded from http://biorxiv.org/on September 18, 2014 



relationships between mutation rates and gene locations (fig. SI and table S2 and S3). 
Comparative analyses showed that patterns of nucleotide substitution rates varied 
substantially among different regions of the genome and among different ethnic 
populations (fig. SI and table S2 and S3). 

These results suggested that some taste receptor variants had different origins 
and emphasized the necessity to narrow the range of target sequences that are to be 
searched on the basis of the ethnic background. Of note, African samples provided a 
huge resource for discovering variant sites, whereas non- African individuals had 
significantly fewer variations (fig. 2C and D and table S2). Furthermore, the 
non-African diversity was largely a subset of the African diversity (fig. 2C and D and 
table S2). 

By comparing alleles among individuals from various ethnic backgrounds, we 
further estimated the extent of differentiation in taste receptor genes before and after the 
divergence from African origins (48) approximately 50,000-75,000 years ago. We used 
a pair- wise proportions test (49), which is used for testing a null hypothesis stating that 
proportions in two populations are the identical. This is referred to as a z-test because 
the statistic is 
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Z = (pr - P2~)/[p~(l " Pi (1 / ni + l/n 2 )] 1,z 

where p" = (pi + p2)/(n i + n 2), and the indices (1,2) refer to the first and second line of 
the table. 

We divided the LoF events of taste receptors into four phylogenies (types A, B, 
C, and D) based on a significant difference (P < 0.05) between populations, suggesting 
that humans might have always been losing their taste receptors, even after population 
differentiation (fig. 2A and E). Moreover, our results revealed that gene mutation flow 
from Africans to non-Africans, and vice versa, had occurred (fig. 2E). This evidence 
suggested that the general theory (50-52) that different evolutionary pressures, such as 
diets, toxins, and climates (energy consumption) have shaped the different 
chemosensory repertoires in mammals might be applicable to modern human 
populations. 

We next used multivariate analysis to investigate country-by-country 
differences in taste LoF variants. To separate confused data sets to make several distinct 
classes, Hierarchical Ward's analysis was used for the matrices of spectra from 14 
ethnic groups to compare their profiles of taste receptor gene variations (fig. 3A, left 
panel; see also fig. S2). Hierarchical Ward's method showed that 14 ethnic populations 
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could be categorized into three general groups: African, Asian, and European-Hispanic. 
Hierarchical median algorithms (fig. S3) and non-Hierarchical clustering algorithms 
(k-means algorithms; fig. 3 A, right panel) also supported these categories. 

Furthermore, we used principal component analysis (PCA) to compare the LoF 
variants of taste receptor genes among various countries (fig. 3B). The three principal 
components (PC 1, 2 and 3) reflected the difference of taste LoF variants in these 
populations. These structure patterns were consistent between the two approaches (fig. 
3A and B). However, when using the country-by-country approach, there were large 
amounts of differences in taste LoF variants among African, European and Hispanic 
populations (fig. 3B). In contrast, taste genetic structures differed only slightly among 
Asian populations (fig. 3B). The taste LoF patterns in Hispanics was genetically closer 
to those in Europeans than the other populations, and some Hispanic (PUR and CLM) 
and African (LWK) groups are similar to Asian populations (fig. 3B). East Asian groups, 
including Japanese and Chinese, had only a few LoF variants of bitter receptor genes 
(table S3). These results illustrated country-by-country genomic signatures for the taste 
LoF variants in humans and sufficient proof that taste receptor genes had evolved 
further in individual colonized areas. 
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Independent origins of mutant alleles 

Previous studies raise the possibility that different ethnic groups often shared 
same taste sensitivity (for example, sweet, salt, bitter and so on) to various compounds 
mediated by taste receptors (53-56). However, phylogenetic relationships among LoF 
taste receptor variants were, for the most part, consistent with the hypothesis that some 
LoF alleles had independently arisen at least twice between the two ethnic groups over 
the course of human evolution (fig. 2A). For example, both rsl50894148 and rs2708381 
caused stop-gain mutations in TAS2R46 at almost the same position, although their 
origins were speculated to be different (rs 150894148: phylogeny type A; rs2708381: 
phylogeny type C; fig. 2A). 

The molecular evolutionary rates of sour and bitter receptor genes seems much 
higher than those of the other receptor genes in any ethnic group (fig. 4). In general sour 
or bitter tastes can be unpleasant, but-because many toxic substances taste sour or 
bitter — they can also be life-saving. My result raised the possibility that modern humans 
lacking sour and/or bitter taste receptors would not seemingly be at a significant 
disadvantage. At present not all of the taste related genes were analyzed, and thus 
confirmation will have to wait until further study is conducted. Characterizing these 
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patterns could facilitate us to unveil the evolutionary pressure acting on modern 
humans. 

Archaic origin of LoF variants in taste receptor genes 

To further address the origin of taste variants, we explored how modern 
humans acquired the highly polymorphic taste receptor genes by analyzing the archaic 
humans Neanderthals (57) and Denisovans (58,59), a likely sister groups to the modern 
human. Because of the vital function, the same sets of taste receptor genes were present 
in both Neanderthals and Denisovans, suggesting that they were not subject to any 
strong multiallelic balancing selection. Virtual genotyping of Denisovan and 
Neanderthal genomes showed great similarities with current humans, though no sign of 
pseudogenization (table S4). Intriguingly, even the most frequent mutation rs 123321 
was not present in Denisovan and Neanderthal genome (table S4). On the basis of these 
findings we concluded that all of the LoF variants in modern human were acquired after 
the divergence of the modern and archaic human and that adaptive introgression of 
population-specific alleles has significantly shaped modern human taste differences (fig. 
5). 
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Genetic diagnoses of LoF variants in taste receptors 

Finally, we performed preliminary genetic testing for LoF variants in using 
author's exome sequences, and examined the effectiveness of the representative list of 
LoF variants described above. The LoF frequency in taste receptors was extremely high 
in the human genome, and this case was no exception: three LoF variants, with two in 
the PKD1L2 gene and one in the TAS2R46 gene (Table 1). Although all mutations were 
heterozygous, we might have a greater risk of vulnerable taste systems with sour and 
bitter sensations. These mutations followed the pattern of the list shown in Figure 2A 
and were unexceptional, which supported the effectiveness of this research. Challenges 
still remain, although our results might yield useful and clinically relevant information 
for the loss of taste receptors. 

Discussion 

In 1931, genetic taste variation was first reported as differences in the ability of 
humans to taste PTC (20). Then many surveys regarding the genetics of taste thresholds 
were conducted to classify PTC tasters and non-tasters roughly according to their taste 
acuity by means of dilutions at which the bitter taste was first detected (21,22). It also 
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has long been known that PTC tasters are nonresponsive to the bitterness in Antidesma 
bunius berries, while PTC non-tasters sense bitterness from these berries (60). Although 
a pile of information about this trait has been accumulated for more than 80 years 
(20-22), its chromosomal location and causative gene were first reported only in 2003 
(13). Over the following ten years, remarkable progress has been made in elucidating 
the genetic variation of taste ability (29-32), but none of the genetic variants explain the 
whole spectrum of interindividual differences of taste sensitivity. 

In this study we first illustrated the widespread occurrence of a recent rapid 

decline in the number of functional taste receptors (pseudogenization). Our report is 

significant because it might raise the hidden possibility that individual taste differences 

may be driven, in part, by LoF variants in taste receptors. On the other hand we also 

think that LoF variants in taste receptors not necessarily mean taste impairment because 

we can consider some scenarios of widespread pseudogenization. Recent study 

suggested that many hTAS2R receptors share a similar chemical ligands and conversely 

one agonist can activate several hTAS2Rs, illustrating that hTAS2Rs cooperatively play 

a common role in toxin avoidance (18,33). Thus some specific taste receptors may have 

been lost in the human as a result of their functional redundancy. An alternative 

possibility is that the structural evolution caused by small amino acid substitutions have 
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made the role of chemical receptors less clear, and taste receptors become unnecessary 
for modern life. Or simply there might be currently no more serious threat to toxins and 
chemical compounds, or multiple receptors may act as backup for gene inactivation. 

Massive loss of taste receptors that we reported is rather unusual, and only a 
few analogous cases have been described. Another pronounced case is olfaction, which 
has also undergone a recent rapid decline in the number of functional genes after species 
differentiation (5,34,35). Olfactory receptors are scattered across the human genome 
and many genes have already been pseudogenized by nucleotide substitutions (5,34,35). 
Olfaction and taste receptors cooperatively contribute to flavor, and thus are speculated 
to have progressed and regressed together during vertebrate evolution (5,34,35). They 
may have evolved by affecting each other with high expansion rates in addition to 
selective pressures due to partial functional redundancy (5,34,35). The high-frequency 
genetic polymorphisms are unique to taste and olfactory receptors in the human genome. 
This co-evolution could comprise a hitherto unexplored aspect of human genotypic 
heterogeneity. Changes in these two sensations provide an important clue with regard to 
future human genetic evolution and may be a major landmark of human evolution. 
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Recent reports provided some examples that non-synonymous mutations cause 
interindividual differences of taste sensitivity (32-35). For instance, Africans have 
higher levels of genetic, geographic and phenotypic diversity at the TAS2R16 and 
TAS2R38 locus (32-35, 61, 62). TAS2R38 P49A, V296I and A262V have clear effects 
on PTC bitterness in humans (32-35). In addition non-synonymous polymorphisms in 
TAS1R1 and TAS1R3 also cause variations of sensitivity to L-glutamate (Umami) 
(29,63). These reports demonstrated that non-synonymous mutations also largely altered 
bitter and umami taste sensitivity among same populations. 

Copy number variants (CNVs) have been reported for TAS2R43 and the -45 
locus (64). CNVs could result in both overrepresentation and absence of expressed 
proteins, and have the potential to exert extreme effects on phenotypes. At present, firm 
conclusions cannot be drawn and more research on genotype-phenotype associations 
should provide insights into the function of each taste receptor. 

My results demonstrated that the LoF variant frequency in taste genes, 
especially sour- and bitter-related receptors, is likely to be extremely high (fig. 2, 4 and 
table S2) although not all of the taste receptors were uncovered. Sour and bitter tastes 
are characteristics of many toxic compounds (65) and have a survival advantage. My 
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result raised the possibility that modern humans lacking sour and/or bitter taste 
receptors would not seemingly be at a significant disadvantage. In ancient times, meals 
that may contain toxic substances were life -threatening for our ancestors (66), but 
because inter-generational or inter-regional exchanges have increased, modern humans 
may have no need to worry about the risks of daily exposures to toxicants. 
Consequently, the human intelligence may have driven LoF variants in taste receptor 
genes. Considering this high mutation rate, "genetically healthy" may be an absolute 
prerequisite for gastronomists and sommeliers. 

I also revealed the country-by-country differences of taste genetic variants. 
African populations had higher levels of genetic diversity at bitter and fat receptors (fig. 
2A and table S3). One possible reason for this might be to increase caloric intake: 
Africans take traditional low-fat, bitter foods, but must eliminate caloric restriction to 
alleviate hunger. In this case LoF variants are beneficial in humans. As shown above 
prevalence of TAS receptors deficiency might be also explained by the hypothesis that 
modern humans may have no need to worry about toxic compounds after acquiring a 
wisdom of living. 
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Recent studies have provided an intriguing example that a loss of sweet 
receptor genes is widespread among carnivorous species and that taste receptors directly 
shape feeding behavior (66-68). In contrast to carnivorous mammals, humans might 
have been losing sour and bitter receptor genes. The recent expansion in number and 
increased survival rates of modern humans are also speculated to be involved in 
high-frequency taste LoF variants. Recent research suggested that 90% variants arose 
within the last 5,000 years (69). 

Many testimonies concerning the taste disorders can be collected in immense 

literatures (70-76). Previously, these disorders could certainly be considered to be 

related to certain conditions such as sequelae of medications (71), radiation therapy 

(71,72), smoking (73), chronic systemic illness (74), and nutritional deficiency (75,76). 

Thus, a strong emphasis has been placed on iatrogenic and secondary taste disorder. In 

addition, various taste disorders are naturally expected to be curable and reversible 

(77-79). My results may raise the hidden possibilities of unknown hereditary taste 

impairment. The mutation frequencies in taste receptor genes could account for the 

number of different modes (70) of specific taste impairment or human tongue -specific 

sensory deficits. It is important to realize that a taste disorder, particularly if it is 

sustained, can have an impact on the quality of a patient's life (80,8). In general loss of 
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taste can lead to nutritional deficiencies due to anorexia, exacerbations of diabetes, 
stroke, and hypertension or even contribute to depression (80,81). These symptoms can 
be particularly serious in older people with chronic illnesses. 

Therefore, we hope that my findings might serve as a basis for clinical 
diagnosis in the near future. Relationships among the taste receptor structures, dietary 
choices, and the associated metabolic pathways remain unclear. Further biochemical 
and cell biological studies are currently underway to uncover the total effect of 
mutations on individual taste sensitivity; they should lead us to a fundamental 
understanding of the relationship between modern human's taste and genetic variants. 

Methods 

Analysis of LoF variants from diverse ethnics 

We analyzed sequence data from 1000 Genomes projects 

(http://browser.1000genomes.org/index.html) and NHLBI (http://www.nhlbi.nih.gov/) 

database. Briefly, these data-sets consisted of high-coverage whole-genome and exome 

sequence data from diverse ethnic groups, respectively. In addition NCBI dbSNP 
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(http://www.nlm.nih.gov/), UCSC genome browser (http://genome.ucsc.edu/) and 
HapMap (http://hapmap.ncbi.nlm.nih.gov/) was used as sequencing platform, which 
were analyzed using an integrated read mapping and variant-calling pipeline to generate 
our initial catalogue of candidate LoF variants of taste receptor genes. MAF (%) was 
calculated using computer simulator. 

Multivariate analysis 

Hierarchical clustering analysis was performed using the R 3.01 statistical software 
together with the Rcmdr package. We employed both Ward's and median algorithms to 
configure the setting for clustering coefficient. Non-Hierarchical clustering algorithms 
was based on k-means approach. To separate confused data sets to make distinct classes, 
principal component analysis (PCA) was also performed on the matrices of spectra from 
14+6 ethnic groups. Two-dimensional score plots and loading profiles of the principal 
components (PC 1 and 2) were applied to visualize the relative contribution of people's 
taste preferences to the clustering of the different spectra. 

Evolutionary scenario of LoF variants in taste receptor genes 
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To estimate how much LoF events in taste receptor genes had occurred before and after 
the divergence from African origins, we compared alleles among individuals of various 
ethnic backgrounds Evidence of phylogeny could be based on significant differences in 
pair-wise comparisons between populations if two groups are significantly different 
(2-sample test for equality of proportions with continuity correction). The standard 
hypothesis test is H>: 7ti = K2 against the alternative (two-sided) H: 7ti ^ K2. The 
pairwise prop test can be used for testing the null that the proportions (probabilities of 
success) in two groups are the same. It is referred to as a z-test because the statistics 
looks like 

Z = (pf - P2WXI - p~) (1 / ni + l/n 2 )] 1/2 

where p " = (p 1 +p 2 ) / (n 1 +n ? ), and indices (1,2) refers to the first and second line of 
the table. In a two-way contingency table where H>: pi = p 2 , this should yield 
comparable results to the ordinary % 2 test. 

Exome sequencing and bioinformatics analysis 

Genomic DNA was extracted from the peripheral blood. Genomic DNA was captured 
on an Agilent SureSelect Human All Exon V5 71 Mb kit (Agilent) following the 
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manufacturer's protocols. DNA was sheared by sonication and adaptors were ligated to 
the resulting fragments. The adaptor-ligated templates were fractionated by agarose gel 
electrophoresis and fragments of the desired size were excised. Extracted DNA was 
amplified by ligation-mediated PCR, purified, and hybridized to the capture array using 
the manufacturer's buffer. Bound genomic DNA was eluted, purified, and amplified by 
ligation-mediated PCR. The resulting fragments were purified and subjected to DNA 
sequencing on the Illumina HiSeq 2000 platform. Captured libraries were sequenced on 
the Illumina genome analyzer as 100-bp paired-end reads, following the manufacturer's 
protocols. Image analysis and base calling was performed by Illumina pipeline versions 
1.4 with default parameters. Sequence reads were mapped to the reference genome 
(hgl9) using the Maq program. Identified variants were annotated based on gene name, 
mutation type and SNP ID. Third-party relatives agreed to a research proposal. The 
risks associated with participation in this research are thought to be minimum because 
the frequency of these variants were high in modern populations. 



Data access 
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Allele information of taste receptors is available from 1000 Genomes projects 
(http://browser.1000genomes.org/index.html) and NHLBI (http://www.nhlbi.nih.gov/) 
database. 
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Figure legends 
Figure 1 

Strategy for identification of deleterious mutations within taste receptor genes 

A flow chart used to identify the sequence variations that have harmful influence on 
human taste senses shows 25 LoF variants within taste receptor genes. Several 
platforms (NCBI dbSNP, UCSC Genes, and Ensemble variation) were used to access 
the validity of variants and examine previous gene annotations. 
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Figure 2 

Loss-of-function variants of taste receptors in major ethnic groups 

(A) Allele frequency of LoF variants of 18 taste genes in European (EA) and African 
(AA) American (NHLBI), and African (AFR), Ad mixed American (AMR), East Asian 
(ASN) and European (EUR) (1000 Genomes). Colored boxes indicate the frequency of 
LoF variants (blue, 0-0.1 %; yellow, 0.1-1 %; green, 1-10 %; red, 10-65 %).Taste genes 
are ordered according to their tastant quality and functions. dbSNP ID and disruption 
type is described. 

(B) Distribution of LoF frequency of each tastant genes. Blue, Green and red lines 
indicate average LoF frequency of taste receptors (2.10 %), olfactory receptors (1.41 %) 
and overall genes of human genome (-0.16 %), respectively. 

(C) Comparison of MAF (%) distribution of taste receptor genes between African and 
European. Vertical and horizontal axis indicate MAF in African and European, 
respectively. Right panel is enlarged view of the red box at the left panel. Blue and red 
dotted circle indicate the population-specific or enriched variants. 

(D) Venn diagram showing the number of shared and unique LoF variants for taste 

receptor genes between European American and African American. 
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(E) Hypothetical origins of different taste between ethnic groups. Divergence of taster 
and nontaster alleles before and after race divergence is indicated as four lines 
(phylogeny type A, B, C, D). 
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Figure 3 

Multivariate analysis of genetic variants across 14+6 ethnic populations 

(A) Hierarchical (dendrogram; ward's algorithm) and Non-Hierarchical (bar graph; 
k-means algorithm) clustering produce three genetic clusters across 14+6 ethnic 
populations (ASW, American's of African Ancestry in SW; CEU, Utah Residents 
(CEPH) with Northern and Western European ancestry ; CHB, Han Chinese in Beijing; 
CHS, Southern Han Chinese; CLM, Colombian from Medellian; FIN, Finnish in 
Finland; GBR, British in England; IBS, Iberian population in Spain; JPT, Japanese in 
Tokyo; LWK, Luhya in Webuye; MXL, Mexican ancestry from Los Angeles; PUR, 
Puerto Rico from Puerto Rica; TSI, Toscani in Italia; YRI, Yorba in Ibadan) (AFR, 
African; ASN, East Asian; AMR, Ad mixed American; EUR, European; EA, European 
American; AA, African American). 

(B) PCA projection of samples taken from a set of 14+6 ethnic groups. Grey-scaled 
spheres correspond to Asian, African, European and Hispanic. 



29 



Downloaded from http://biorxiv.org/on September 18, 2014 

Figure 4 

Molecular evolutionary rate of five taste senses in modern humans 

Average LoF mutation frequency of taste receptor genes of five senses (sour, bitter, 
umami, sweet and salt) are shown. High and low groups are presented according to the 
molecular evolutionary rate. 
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Figure 5 

Evolutionary history of loss-of-function variants of taste receptors 

Schematic tree shows LoF events of taste receptor genes in modern and archaic humans. 
Red arrowheads show mutation events. Dashed lines indicate the major introgression 
between ethnic groups. 
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Table 1 

Genetic diagnosis of author's taste receptors 



Gene name 


Transcript ID 


Chromosome 


Position 


dbSNP 


Mutation type 


Zygosity 


ref peptide 


Coverage 


PKD1L2 


NM_052892 


chrl6 


81183325 


rs59980974 


Stop Gain 


Hetero 


K/* 


28 






chrl6 


81199544 


rs 12925771 


Stop Gain 


Hetero 


R/* 


22 


TAS2R46 


NM_176887 


chr!2 


11214145 


rs2708381 


Stop Gain 


Hetero 


W/* 


57 
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rsl 508941 48 stop-gained 
rs2708381 stop-gained 
rsl 38984388 stop-gained 










































































































































































Sodium ENaCd 


rsl 48588796 stop-gained 














CD36 rsl 14881528 5'splice 
p . rs3211893 5'splice 
rat rs3211938 stop-gained 

rs56381858 stop-gained 



















































0.1-1% 



B C 
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Molecular evolutionary rate 

Fast Slow 



Sour (14.7%) 
Bitter (1 .8%) 



Umami (0.17%) 
Sweet (0.02%) 
Salty (0.05% ) 
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Mutation 




no Mutation 



